Entanglement signatures of quantum Hall phase transitions 
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We study quantum phase transitions involving fractional quantum Hall states, using numerical 
calculations of entanglements and related quantities. We tune finite-size wavefunctions on spherical 
geometries, by varying the interaction potential away from the Coulomb interaction. We uncover 
signatures of quantum phase transitions contained in the scaling behavior of the entropy of entan- 
glement between two parts of the sphere. In addition to the entanglement entropy, we show that 
signatures of quantum phase transitions also appear in other aspects of the reduced density matrix 
of one part of the sphere. 



I. INTRODUCTION 

An exciting interdisciplinary development in recent 
years has been the description of condensed matter 
phases using entanglement measures borrowed from the 
field of quantum information theory One very re- 
cent example is the characterization of topological order 
in fractional quantum Hall states 0, H, |j| using entan- 
glement. Topologically ordered ground states are char- 
acterized by fractionalized excitations, degeneracies on 
higher-genus surfaces, and an energy gap in the excita- 
tion spectrum above the ground state For exam- 
ple, on a genus-5 surface, the Laughlin state at filling 
v = l/m has ground states. Although a number of 
spin models theoretically possess topological order, the 
only confirmed experimental realizations of topological 
order are fractional quantum Hall (FQH) states of two- 
dimensional (2D) electrons in a magnetic field. In this 
Article we focus on this most realistic class of topolog- 
ically ordered states. FQH states have been the object 
of further intense attention because of the possibility of 
cold-atom realizations [6*1 , and more recently due to quan- 
tum computation proposals based on their topological 
properties Q- 

In considering entanglement properties, we will focus 
on bipartite entanglement between two parts (A and B) 
of the system. This is characterized by the reduced den- 
sity matrix pA — trs p oi subsystem A, obtained by trac- 
ing out all the B degrees of freedom. While it is instruc- 
tive to study the upper part of the eigenvalue spectrum of 
Pa, the so-called entanglement spectrum, it is also often 
convenient to extract a single number from this spectrum. 
For the latter purpose we will use the entanglement en- 
tropy, Sa ^ -tr [pA In pa] ■ 

To study topological order in FQH states, Refs. [1, Q 
exploited the novel concept of topological entanglement 
entropy which appears in the entanglement entropy 

Sa between a block A and the rest (B) of the system. 
On the other hand, Ref. Q studies the spectrum of pA 
and extracts information relevant to the topological order 



from one end of the spectrum. 

A quantity characterizing a phase should also provide 
signatures of quantum phase transitions leading into or 
away from that phase. This idea is already fruitful in ID 
physics, where the (asymptotic) block entanglement en- 
tropy can distinguish between gapless and gapped phases. 
This allows one to pinpoint the phase transition between 
gapped and gapless phase, in particular in DMRG calcu- 
lations where the block entanglement entropy is readily 
available. For the case of quantum phase transitions in- 
volving topologically ordered states, entanglement stud- 
ies have been exploited for the simpler 'toy' case of Kitaev 
models with order-destroying additional terms ITl| . 
In this Article, we will examine the utility of entangle- 
ment calculations for studying more 'realistic' topologi- 
cal phase transitions, namely, transitions involving FQH 
states. FQH states being the only experimental exam- 
ples of topological order, this is an important step toward 
developing entanglement-based tools that are useful for 
experimentally relevant situations. 

We study phase transitions between FQH and non- 
FQH states driven by a change in the interaction poten- 
tial. Specifically, using the Haldane pseudopotential de- 
scription of interactions projected to specific Landau lev- 
els [12j , we vary the first pseudopotential Vi away from its 
Coulomb potential value. The transition we mainly focus 
on involves the best-known FQH state, for fermions at 
filling ly = 1/3. With a Coulomb interaction between the 
fermions, the ground state is known to be an FQH state 
topologically equivalent to the Laughlin state [iSj] . If the 
the first pseudopotential is reduced enough, the Laughlin 
state ceases to be a good description of the ground state. 
There is thus a phase transition as a function of Vi [l3| ■ 
We will present entanglement calculations which probe 
this phase transition. We also present results for filling 
fraction v — 5/2, where the possibility of a more intricate 
quantum Hall state (the Moore-Read state [l!]), provides 
a more challenging situation. 

In Sectionini we summarize necessary background ma- 
terial on the topological entanglement entropy 7 and on 
FQH wavefunctions and transitions, and develop the seal- 
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FIG. 1: (Color online) Possible behaviors of so, defined in 
Eq. IT}, within a topologically ordered phase {X < Xc) and 
after a phase transition into a different phase {X > Xc). 



ing concepts necessary for treating phase transitions. In 
the subsequent sections, we use three different aspects of 
entanglement to probe these phase transitions. First, we 
consider the entanglement entropy of a block with the 
rest of the system, and track phase transitions through 
the behavior of this quantity as a function of block size 
and system size. Since our calculations are based on 
finite-size wavefunctions, there are limiting procedures 
involved which can be performed in different orders. Re- 
sults obtained using different limit orders are discussed 
in Sections inilllV] and El 

Second, we consider the top part of the reduced density 
matrix spectrum. This analysis is based on the identifi- 
cation of features of the entanglement spectrum in terms 
of topological order and related edge physics 3- We 
show in section IVII how the entanglement spectrum is 
affected when the system is driven through a quantum 
phase transition. 

Finally, in section IVlIl we use the concept of majoriza- 
tion, which is based on comparing complete reduced den- 
sity matrix spectra for two wavefunctions. Majorization 
relations between reduced density matrices obtained from 
condensed matter wavefunctions has been the subject of 
intriguing recent studies [H, 113. While the full impli- 
cations are not yet clear, this work adds to the growing 
understanding of majorization in condensed-matter sys- 
tems. 



II. ENTANGLEMENT AND FQH STATES 

A. Block entanglement entropy and topological 
phase transitions 

For 2D topologically ordered systems, an important 
recent result [1, Q relates the block entanglement entropy 
Sa to the topological order. The dependence of the block 
entanglement entropy Sa on the length L of its boundary 
is asymptotically linear, in accord with the "area law" . 



In addition, Refs. fs!, [§] have found that this dependence 
also has a topological sub-leading term: 



Sa{L) 



aL 



7 



The sub-leading term is called the topological entan- 
glement entropy (abbreviated topological entropy) and 
is a constant for a given topologically ordered phase, 
7 = InP. Here T) is the total quantum dimension charac- 
terizing the topological field theory describing the phase. 

The concept of topological entropy is an important new 
development, because the usual definition of topological 
order is often unwieldy to use in practice for theoretical 
or numerical studies of such order. The entanglement- 
based characterization provides a new route for iden- 
tifying topological order and, by extension, topological 
phase transitions. After the quantity (7) was introduced 
for topologically ordered states in general, it has been 
explored in several specific contexts, including quantum 
Hall states 0, H, [l8l[l9l[23l , quantum dimer models [2l[ , 
and Kitaev models jlfl [nf, [13, [11, [11] . 

We will employ considerations of both the leading lin- 
earity and the sub-leading invariant term to characterize 
topological phase transitions. Let us consider, very gen- 
erally, quantum phase transitions between a topologically 
ordered phase and another phase. For the block entan- 
glement entropy, let us imagine that we have determined 
the asymptotic relationship 
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where L is the boundary of the block. Note that this is 
not always possible; in some 2D phases the leading term 
might not be purely hnear [H, [ll, [13, [11] . Note also that 
the above behavior does not necessarily imply topological 
order; Ref. [1^ gives an example of a non-topological 
state following such a relationship with nonzero sq. 

If we are in a topologically ordered phase, the negative 
intercept sp will by definition be equal to the topological 
entropy, 7 = In I?. Figure [T] shows some possibilities of 
what can happen to sq as the 2D system is driven across 
a quantum phase transition away from the topologically 
ordered state, by varying a parameter X across the criti- 
cal value Xc- In the parameter range X < Xc where the 
system is in the topologically ordered phase, sq is fixed 
at a positive plateau (sq = 7). 

Case A shows a transition into another topologically 
ordered state with a different quantum dimension; sq 
jumps to another constant value 7'. The other figures 
show transitions to non-topological phases. Case B shows 
a transition to a gapped state which is not topologically 
ordered - the intercept so drops to zero. (Ref. 10] treats 
an example.) Case C shows a transition into a non- 
topological phase, in which sq is nonzero but not con- 
stant. Finally, Case D shows a transition into a state 
where the leading term in the asymptotic behavior of 
Sa{L) is not linear, so that sq is undefined. 



3 



B. FQH wavefunctions on spheres 

For finite-size studies of fractional quantum Hall 
physics, spherical and toroidal geometries are particu- 
larly popular when one wants to avoid complications due 
to the edge. We set notation by providing a rapid review 
of FQH wavefunctions on spheres [H, [3] ■ We study N 
electrons on a sphere of radius i?, subjected to the mag- 
netic field B provided by a magnetic monopole at the 
center of the sphere. The Dirac quantization condition 
requires the number of fluxes — At:B?B /(j)^ to be in- 
tegral, where ipQ = h/eis the quantum of flux. If we mea- 
sure length in units of the magnetic length I b = yh/eB, 
then the quantization condition is i? = yJ~N\fj2. The 
1 Landau-level orbitals are labeled either Z = to 



v = 1/3 



v = 5/2 





Nfj, or Lz — —L to +L. The wavefunctions in spherical 
coordinates are given by 



\MO,v)\'^{cos- 
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There is no dependence on the azimuthal angle; the or- 
bitals are each localized around a "circle of latitude", 
with the / = orbital localized near the "north pole." 
Basis states for the A^-electron wavefunctions are ex- 
pressed in terms of the occupancies of these orbitals. 

To define entanglements between two parts of the 
sphere, one has to first partition the sphere into A and B. 
We define block A to be the first Ia orbitals, extending 
spatially from the north pole out to some latitude, i.e., 
including orbitals I = through I = Ia — I- 

While neighboring orbitals do overlap, it is natural to 
think of the location where their amplitudes are equal, 
I ^f, (61) 1 2 = |^',+i(6')|2, as the "boundary" between the 
orbital I and the orbital {I + 1). This happens at angle 



9 oi 
tan^ -7- = 



l + l 



I 



The boundary between partitions A and B is thus a circle 
of latitude at polar angle di^-i- The boundary length is 



CiAN) 



2ttR sin 6*; ._i 
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In finite-size spherical geometries, FQH states do not 
appear exactly at the filling fraction v. Instead, the num- 
ber of orbitals (iV^ -I- 1) is related to the number of par- 
ticles N by: 



S 



where the shift S is an integer determined by the FQH 
state [29l|. Although S is insignificant in the thermody- 
namic limit, it can create an "aliasing" problem in nu- 
merical studies: different FQH states can compete at the 
same finite values of N and N^. 
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FIG. 2: (Color online) Overlaps of ground state wavefunc- 
tions at interactions relevant to fillings 1/3 and 5/2, respec- 
tively with Laughlin and Moore- Read wavefunctions. A large- 
and a small- case are shown for each fraction. 



C. Exact diagonalization 

Our analysis is based on numerical exact diagonaliza- 
tion of Coulomb-like Hamiltonians projected onto appro- 
priate Landau levels (lowest Landau level for u — l/Z and 
the second Landau level for v = 5/2). We employ the 
spherical geometry described above. For this work, we 
have calculated = 1/3 wavefunctions up to 14 fermions 
and 1/ = 5/2 wavefunctions up to 20 fermions. To give 
an idea of the scale of these calculations, we note that for 
the = 20 wavefunctions at v = 5/2 the Hilbert space 
dimension is 193498854, while for the A^ = 14 wavefunc- 
tions dX V = 1/3, it is 129609224. Dimensions are given 
here after reduction using the discrete symmetry under 
global Lz flip. 

We used the Lanczos algorithm to compute the sys- 
tem ground state. Calculations were performed on a PC 
cluster with 24 cores (AMD Opteron 265) and 48 Gbytes 
of memory. A Lanczos iteration typically requires up to 
17h of cpu time for the largest Hilbert spaces. In our 
study of phase transitions, calculations are particularly 
time-intensive when the system is in or near a gapless 
phase, in which case the Lanczos procedure requires a 
larger number of iterations to converge. 



D. Phase transitions 

For FQH systems, the natural interaction parameter to 
vary is one of the Haldane pseudopotentials Vm obtained 
by decomposing the interaction potential into channels 
specified by the relative angular momentum ni of the in- 
teracting particles [l3, . We use interaction potentials 
whose m = 1 pseudopotential Vi is changed, while the 
other Vm are fixed at the Coulomb value for that Landau 
level. We calculate and present wavefunction properties 
as a function of 5Vi cx [Vi — V™"'). Tuning of Vi loosely 
represents variable aspects of FQH experiments, such as 
the thickness of the quantum well where the 2D electron 
gas resides. 

For the v = 1/3 case, the exact Laughlin state is 
obtained for 5Vi — > -l-oo (short-range interaction), and 
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continues to be a good description of the state for the 
Coulomb potential, 6Vi = 0. However, for negative 
SVi cx {Vi — V7^°"'), there is a quantum phase transi- 
tion at some dVi — 6Vic < into a non-FQH incom- 
pressible state The overlap plots of Fig. [2] (left 
panel) shows the transition to be somewhere between 
5Vi = -0.065 and SVi = -0.11. The non-FQH state 
for 5Vi < 6Vic presumably has charge-density-wave or 
Wigner-crystalline order [1J|; the details are not impor- 
tant for our purposes. 

For — 5/2, the overlaps (figured]) suggest that the 
Moore-Read state is stable in some window of slightly 
positive SVi (around SVi ~ 0.03), and that there are 
transitions on either side of this phase to non-FQH phases 
[sol . [Slj . The non-FQH phases are possibly a striped 
charge-density-wave phase on the left of the Moore-Read 
region and a composite fermi sea on the right [Slj . 
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The definition of the topological entanglement entropy 
7, or the quantity sq in equation [Tl involves two extrap- 
olations: (1) to the thermodynamic limit, N oo, and 
(2) to the asymptotic limit of the size of the A block, 
i.e. Ia ^ ^- This double scaling limit can be approached 
via different possible extrapolation paths in the {N, I a) 
space. To describe the different extrapolation methods, 
we first rewrite the relation ^ for finite N and Ia- 



FIG. 3: (Color online) Top four: Si^{N) versus 1/iV for 
several SVi values, displaying how the smooth behavior of the 
these curves disappear. Bottom left: quality of fit (x^) to a 
simple function, Co -|-ci/A^^. Bottom right: Attempt to extract 
So. The quantity is less and less meaningful for more negative 
SVi (discussion in text). The topological value ilnS « 0.55 
is shown by dashed horizontal line. 



F. The entanglement spectrum 



Si AN) = s.CiAN) - so{N) (2) 

We have ignored ^^i-dependence of sq and both N- and 
Zyi-dependence of si. Our experience with numerical 
wavefunctions indicate that these dependences are gen- 
erally weak, at least for model FQH wavefunctions. Note 
that, even for topologically ordered model FQH states, 
the block entanglements Si ^ are not linear in \/Ta for fi- 
nite N [c.f. Fig. 1 in Ref. [2]). The reason is that ^/Ta 
is proportional to the square root of the area of the A 
region, which is not proportional to the circumference 
C/^ {N) in the curved geometry of the sphere. 

In Sec. IIIH we consider the method used in Refs. 0, 
namely, performing the N ^ oo extrapolation first 
for each I a, and then using the resulting iS';^(oo) versus 
C/^(oo) ^ \/Ta to extract so(oo) or 7. 

In Sec. IIVI we consider the extrapolation procedure in 
reverse order, namely, first extracting so{N) from finite- 
size Si^ versus C;^ dependences, and then taking the 
iV — > CX) limit. 

In Sec. |Vl we take the block A to be half the system 
(which is a natural choice in a spherical geometry) , I a = 
l*^ = l(N^ + 1). Examining the Si-^N) versus CrjN) 
may be regarded as taking the two limits simultaneously. 



Clearly, the complete spectrum of the reduced den- 
sity matrix pA of a subsystem A contains more informa- 
tion than any one number (such as the entropy of entan- 
glement Sa) extracted from this spectrum. Extraction 
of additional information from the complete spectrum 
has been reported in several other condensed-matter con- 
texts. 

Ref. Q has empirically shown that, for FQH wavefunc- 
tions on a sphere, the spectrum of the reduced density 
matrix of one hemisphere can be related to the confor- 
mal field theory (CFT) describing the ID edge of the 
FQH state. We exploit this notion in Sec. IVII to show 
how quantum phase transitions appear in the so-called 
entanglement spectrum. 

Another way of exploiting the complete spectrum of 
Pa is via majorization comparisons, which we use in Sec. 

ivni 

III. EXTRAPOLATION FOR EACH BLOCK 
SIZE 

We attempt to employ the extrapolation method of 
Refs. 0, 3, first extrapolating N ^ 00 ior fixed values 
of Ia- From the Si^{N) versus 1/N plots in figure[3](top 
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4 panels), we note that the extrapolation starts to lose 
meaning as one reduces Vi beyond the presumed transi- 
tion. This is made quantitative by estimating goodness- 
of-fit of the Si^{N) versus 1/N data to simple func- 
tions. The plotted estimates are obtained from trying 
Sij^{N) = cq+ci/N'^; any other reasonable function gives 
similar results. This versus SVi curve can be regarded 
as one entanglement signature of the phase transition. 

The jump in gets sharper for smaller I a- In partic- 
ular, for = 1 it is distinctly localized at 5Vi ~ —0.1. 
This appears to provide a sharp estimate for the location 
of the transition. One tempting interpretation is that, 
since the system is large compared to a single orbital, 
Sij^=i{N) for moderate N is already a good indicator of 
the thermodynamic limit; hence the sharp jump. 

One could still estimate 5;^ (TV oo) extrapolates, 
disregarding the "noise" in the Sij^{N) versus dat- 
apoints. The resulting S'i^(oo) points, plotted against 
^/Ta, gives a line whose intercept is by definition —sq. 
The estimates of sq, thus obtained, move away from the 
expected value of 7 = In -^/S ^ 0.55 around the transition 
point. However, it is important to note that in the region 
6Vi < —0.06, the estimates have relatively little mean- 
ing. First, there is the increasing scatter in the SlJ^{N) 
versus 1/7V data, which makes the S';^(cxd) values unreli- 
able. Second, we find that the S';^(oo) versus ^/Ta curves 
are more and more curved as one increases —SVi further 
away from the Laughlin state. Given the uncertainties, 
we do not attempt to estimate error bars for sp- The 
point here is to note how sq loses meaning, in the large- 
(—SVi) region for which we have used shaded symbols in 
figure [3l Referring back to figure [H the situation we have 
should be regarded more like case D and not like case C. 



IV. DEPENDENCE ON CIRCUMFERENCE 
AND AREA LAW 

We now consider extrapolations in reverse order, i.e., 
first find a negative intercept so{N) for each N by looking 
at the large-Z^ behavior at that N, and only afterwards 
consider the A^-dependence. For topologically ordered 
model FQH states, the block entanglements 5;^ (N) have 
good linear dependences on the circumference Ci^ (N) oc 
\/Ia{N^ + 1 - Ia), for any fixed A^. 

Here we are of course interested in realistic wavefunc- 
tions, i.e., the ground states of Coulomb-like potentials. 
In figure HI Si^ versus C/^ plots are shown for several 
SVi, for = 1/3 states of 12 particles. It is interest- 
ing to note that the plots remain smooth as —SVi is in- 
creased past the transition; however they acquire curva- 
ture. As the dependence deviates from the linearity of 
model wavcf unctions, extracting sq loses meaning. We do 
not attempt estimating sq. Instead, the deviation from 
linearity is shown through measures. The rise in 
again represents the quantum phase transition. 

Note that this deviation from linearity does not nec- 
essarily constitute a violation of the area law, which is 
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FIG. 4: (Color online) Top panels: entanglement entropy 
plotted against circumference of A block; v = 1/3, A*' — 12. 
The linearity is progressively destroyed as one reduces the 
pseudopotential Vi . Lower panel: for linear fits to the S;^ 
versus circumference data. 



a statement about the thermodynamic limit. The devi- 
ation could instead indicate a stronger ^-dependence of 
the quantity si in Eq. 



ENTANGLEMENT BETWEEN EQUAL 
HEMISPHERES 



We now consider entanglement between equal 
hemisphere-shaped A and B blocks, i.e., we use I a = 
l*^ = {N^ + l)/2. According to Eq. ©, SrJN) should 
depend linearly on the circumference C;^ (N) = 2-kR cx 
^J Ntf,. For the exact Laughlin state, an example is shown 
in Fig. 1 (inset) of ReL [2]. Here, we explore the fate of 
this linearity as a function of 5Vi . 

The top row in figure [5] plots Si*^ against y/N^j, for 
several Vi values at u = 1/3. The points acquire scat- 
ter as —5Vi is increased, until it becomes meaningless to 
extract the intercept; the situation is similar to Section 
mil (figure [3|) . This process is indicated by the curve 
(center panel) signifying the goodness of the linear fit. 

The 1/ = 5/2 case (figure [5] bottom panel) similarly 
shows that the scatter is low (so that a linear Si* versus 
fit is meaningful) only in the window of 5Vi where 
the Moore- Read state describes the physics of the ground 
state. 

A technical note: for even 7V^, the number of orbitals 
is odd and it is impossible to divide the sphere precisely 
into halves. In this case we keep 1*a = orbitals in 

block A, which nearly divides the sphere into halves, and 
use values Si*^=m^/2 for the entropy and C/^=7v^/2(A^) for 
the circumference. 
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FIG. 5: (Color online) Top row: entanglement entropies for 
hemisphere partitioning, plotted against sphere radius; i/ — 
1/3. Center row shows quality of linear fit as function of 
pseudopotential SVi for u = 1/3. Bottom row shows for 
ly = 5/2. 



VI. ENTANGLEMENT SPECTRUM AND 
ENTANGLEMENT GAP 
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FIG. 6: (Color online) Top panels: entanglement spectrum, 
1/ — 1/3, A'' = 12, block A containing Ia = 17 orbitals and 
Na = 6 fermions. Main plot: ground state for SVi = 0.04. 
Ellipse indicates the most prominent "conformal" part of the 
spectrum. Arrow indicates the "entanglement gap" So be- 
tween CFT and non-CFT parts of the spectrum. Inset shows 
exact Laughlin state, which has no higher-lying non-CFT 
part. Lower panel: Empty dots show two lowest levels at 
Z/i'''' = 54, plotted against SVi. Filled squares show "entan- 
glement gap", the difference of the two lowest levels. 



Following Ref. we introduce the "entanglement 
spectrum" ^ as Ai = exp(— ^i), where are eigenvalues 
of the reduced density matrix pA of one hemisphere. The 
eigenvalues can be classified by the number of fermions 
Na in the A block, and also by the total "angular mo- 
mentum" Li"^-* of the A block. It was argued [3] that 
the low-lying spectrum of the reduced density ma- 
trix for fixed Na, plotted as a function of Li^\ should 
display a structure reflecting the conformal field theory 
(CFT) describing the edge physics. In figurelHlthis "CFT 
spectrum" is marked with an ellipse. For interactions at 
which the FQH state provides a good description of the 
physics, the CFT spectrum is well-separated by a gap 
from a higher "non-CFT" part of the spectrum. 

As in Ref. [3, we denote the gap between the lowest 

(A) (A) 

two , at the Lz value where the highest-Lz member 
of the CFT spectrum occurs, as 6o- In figureEl this is the 
gap between the lowest two states at ii^' — 54 (marked 
by arrow). We study what happens to the spectrum as 
we tune the interaction away from the FQH state across 
a quantum phase transition. We quantify the change of 
the spectrum in terms of the quantity Sq, defined above. 



(The quantities 61,2 defined in Ref. [3], the gaps at other 

(A) 

Lz values, are expected to have similar dependence on 
SVi.) 

In figure [S] (lower panel), we plot Sq as a function of the 
pseudopotential SVi for the 1/ = 1/3 case. This clearly 
shows a dramatic decrease of the "entanglement gap" 
around the region of the phase transition. The two levels 
in question are also individually plotted with open dots; 
there is a level crossing around SVi ^ —0.1. We note 
that for values of SVi < —0.1 the CFT-like structure 
of the entanglement spectrum is lost so it is no longer 
meaningful to think of Sq as the gap between CFT and 
non-CFT energy levels. A similar picture is observed for 
Moore- Read wavefunctions ^ . 



VII. MAJORIZATION 

The concept of majorization involves comparison of 
two complete spectra. In the context of condensed- 
matter applications, it generally involves the comparison 
of two reduced density matrix spectra corresponding to 
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FIG. 7: (Color online.) Left panels: Entanglement entropy 
and majorization plotted against SVi for — 12, i/ = 1/3. 
Right panels: same for A'' = 18, v — 5/2. Majorization at 
some SVi is taken to be 1 if the corresponding pA spectrum 
majorizes a spectrum at an adjacent value of SVi; otherwise 
it is set to 0. 



the ground states of two Hamiltonians with shghtly dif- 
ferent parameters [H, [13, HI] . 

To define majorizaton, we consider two sets of n real 
elements {A^} and {/Lt^} sorted in decreasing order and 
satisfying J2i = Si Mi — 1- One says that the set {A^} 
majorizes the set {^i} if 



Vk € {1, 



This relationship is often expressed as A /i. Obviously, 
if A >- /i then S{X) < S{fi), where S{x) = — Xi Inxi is 
the von Neumann entropy. 

In ID quantum systems, the spectrum of the reduced 
density matrix for a spatial block has been argued [H, 
[TtI . [ssj to become more majorized as one moves along a 
renormalization group (RG) flow, away from an RG fixed 
point. As of now, there are no established general results 
for 2D quantum states. 

In figure [71 we examine majorization of reduced den- 
sity matrices between ground states at different values of 
the pseudopotential Vi . We show results for the v — 1/3 
system with iV = 12 particles, and the = 5/2 system 
with TV = 18 particles. For the i' — 1/3 system, we find 
that the eigenspectrum is continuously majorized as SVi 
decreases down to the value SVi = —0.10, i.e., down to 
the phase transition region. In this region, the eigenspec- 
trum of Pa for each SVi majorizes the pA eigenspectrum 
at a smaller (more negative) 6Vi. For SVi < —0.1, the pA 
spectra are not majorized. This result is robust for differ- 



ent sizes {Ia) of the A block. A similar picture emerges 
for the V = 5/2 system; majorization occurs in a region 
roughly corresponding to where the ground state has the 
structure of the Moore-Read state. However, the effect 
is more fragile, e.g., the extent of dVi values where the 
majorization is found, depends on the partition size {I a) 
used. Note for v = 1/3 filling that Si^ has a kink near the 
transition point. No such feature is seen for the v = 5/2 
case, neither in Si^{Vi) not in its derivative. 

To summarize, we have demonstrated for v — 1/3 that 
the region where majorization occurs coincides dramati- 
cally with the region where the system is in an FQH state; 
we no longer find majorization away from this phase. The 
situation is similar but less clear ior v — 5/2. A full un- 
derstanding is lacking at the moment, but several intrigu- 
ing speculations present themselves. Most obviously, it 
is tempting to think of majorization being an indicator 
of quantum phase transitions. Second, since the reduced 
density matrix of a block of the sphere contains informa- 
tion about the physics of the quantum Hall edge Q , it is 
possible that our majorization results can be interpreted 
in terms of the evolution of the edge as a function of 5Vi . 



VIII. CONCLUSIONS 

We have presented a numerical study of how quan- 
tum phase transitions involving fractional quantum Hall 
states are manifested in entanglement measures and re- 
lated quantities. We used three extrapolation methods 
to examine the double-scaling limit where the block en- 
tanglement intercept (sq in equation [T]) is defined. We 
showed that the breakdown of these extrapolation proce- 
dures signals the quantum phase transitions away from 
the topologically ordered incompressible FQH states. In 
addition, the entanglement spectrum was used in more 
detail in two different ways to characterize the quantum 
phase transitions, exploiting recently-developed concepts 
of CFT spectrum and majorization. 

Our work opens up several new open questions and 
directions. 

Since we have explored topological phase transitions on 
spherical geometries only, it would be instructive to look 
for analogous signatures on a toroidal geometry. Entan- 
glement scaling behaviors are known in much less detail 
for FQH states on the torus, and further investigations 
are clearly needed. 

Second, our results on majorization invite a thor- 
ough investigation of this concept, in general for two- 
dimensional systems and in particular both for 2D phase 
transitions and for topologically ordered states. It would 
appear that the knowledge necessary for putting our find- 
ings in context does not yet exist. The study of majoriza- 
tion, and the upper part of the reduced density matrix 
spectrum, also raises the question of other signatures of 
topological phase transitions one might yet extract from 
the full Pa spectrum. 

Finally, for topologically ordered states, it is promising 
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to explore our three extrapolation methods, for extract- 
ing the topological quantity 7. In previous work with 
FQH states on spherical geometries [3, |3| , the focus has 
been on the first extrapolation procedure (section |TTT|. 
If the topological entanglement entropy is to be used to 
identify intricate FQH states and their CFT's, improved 
methods of calculating 7 are vital. Alternate extrapola- 
tion methods is one direction that needs to be explored 
in this regard. 
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